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Abstract 

In this paper the study of "7 + Jet" events is continued, aimed at jet energy scale setting and hadron 
calorimeter calibration at LHC energies. The event number distribution over P 4 7 and rp is presented. The 
features of "7 + Jet" events in the barrel region of CMS detector (|?7| < 1.4) are exposed. P t structure of 
the region in the 77 — <fi space inside and beyond jet is also shown. 



1. INTRODUCTION 



In this paper we continue the study of "7 + Jet" process, started in [1] aimed at jet energy 
scale setting and hadron calorimeter calibration at LHC energies. 

This article is organized in the following way. In Section 2 we shall estimate the 
event number distribution dependence on a value of "back-to-back" <f»/~ j et ) angle between 

-» 7 -> Jet 79 R 

the direct photon P t and jet P t and on the P t of the initial state radiation (ISR), i.e. P t 
value. PC and rf dependence of rates is also studied in this section. Estimation of "7 + 
Jet" event rates for Barrel (HB), Endcap (HE) and Forward (HF) parts of hadron calorimeter 
(HCAL) is performed in the last subsection of Section 2. In Section 3 we consider the "7 + 
Jet" event features in HB region. In particular, we study the Ptcur parameter influence 
on the selection of events with a good photon and jet P t balance and also the problem of 
selection of events with suppressed ISR activity. The next subsections are devoted to the P t 
structure of jets in 77 — <fi space and to the particles P t distribution inside and beyond the jet. 



2. EVENT RATES FOR DIFFERENT Pp AND rf INTERVALS. 
2.1 Event number distribution dependence on "back-to-back" 0( 7 je t) angle and on 
P t ISR values. 

The definitions of observable physical variables introduced in Section 3.2 of [1] allow us to 
make the first step in the study of a possible way to select the events with a good Pt 1 and p t Jet 
balance. Here we shall study the spectrum of Pt 56 variable (that is proportional to Pt 
up to the value of intrinsic parton transverse momentum inside a proton) in signal events 
with direct photons and this spectrum dependence on some of experimentally measurable 
variables. For this reason the results of 10 7 "7 + Jet" events generation for every of four 
Pt 1 intervals (40-50, 100-120, 200-240, 300-360 GeV/c) by use of PYTHIA with account 
of 2 QCD subprocesses qq — > g + 7 and q g —> q + 7 were analyzed. For this analysis 
the observables and the "Selection 1" cuts (1-7), defined in Sections 3.2 of [1], with the 
following cut parameters: 

PtwT = 5GeV/c; e^or = 7%; < 15°; P t %l% = 20 GeV/c. (1) 

were used. The value of the luminosity was taken to be L = 10 33 cm 2 sec _1 (low LHC 
luminosity). 

In Tables 1, 2 and 4, 5 we present P t 56 spectra for two most illustrative cases of 
Pt 1 intervals: 40 < PP < 50 GeV/c (Table 1 and 2) and 200 < P t 7 < 240 GeV/c 
(Table 4 and 5). The distributions of the number of events for integrated luminosity Li n t — 
3 in different intervals of Pt 56 {(kr) was taken to be fixed at PYTHIA default value, i.e. 
(kr) = 0.44 GeV/c) and for different "back-to-back" angle intervals 4>(-y,jet) = 180° ± A<fi 
(A(f> = 15°, 10° and 5° as well as without any restriction on A(f>, i.e. for the whole <f> interval 
Acj) = 180°) [] are given there. Tables 1 and 2 correspond to the case p t clust < 20 GeV/c 
and serve as an illustration since it is rather a weak cut condition, while Tables 4 and 5 
correspond to a more restrictive selection cut value Ptcur = 5 GeV/c. 

Firstly from the last summary lines of Tables 1-4 we can make a general conclu- 
sion about the A</>-dependence of the events spectrum. Thus, in the case of weak restriction 

'The value A0 = 5° coincides with one CMS HCAL tower size in the </>-plane. and, thus, defines the precision 
of experimental registration of the jet Pt vector. 
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pdust < 20 Ge y/ C we can see from Table i that for the interval of 40 < P t 7 < 50 GeV/c 
there are about 66% of events concentrated in the A</> < 15° range, while 32% of events are 
in the A(j> < 5° range. At the same time the analogous summary line of Table 2 shows us 
that at higher values 200 < Pt 1 < 240 GeV/c the events spectrum moves to the small A<fi 
region: more than 99% of events have Acf> < 15° and 79% of them have Acf> < 5°. 

A tendency of very rapid concentration of the signal "7 + Jet' event distributions 
with the growth of P t 7 in a rather narrow back-to-back angle interval A<fi < 15° becomes 
more distinctive with a more restrictive cut Ptcur = 5 GeV/c (see Tables 4 and 5). From 
the last summary line of the Tables 3 we see that in the case of 40 < P t 7 < 50 GeV/c more 
than 97% of the events have A<fi < 15°, while 68% of them are in the A<f> < 5° range. 
At 200 < Pt 1 < 240 GeV/c (see Table 5) more than 99% of the events that undergo this 
cut have A(f> < 5°. It means that imposing Ptcur = 5 GeV/c and suppressing cluster or 
mini-jet activity, we choose the events with a clean "back-to-back" (within 15°) topology of 
7 and jet orientation. 

So, one can conclude that PYTHIA simulation predicts that at LHC energies most of 
the "7 + Jet" events after imposing Ptcur — 20 GeV/c may have the vectors Pt' and Pt 
being back-to-back within A(j> < 15°. The cut Ptcur — 5 GeV/c improves this tendency [| 

It is worth mentioning that this picture reflects the predictions of one of the generators 
which are based on the approximate LO values for the cross section. It may change if the next- 
to-leading order or soft physics effects are included and also after more real experimental 
data would be collected for comparison with a theory. 

The other lines of the same Tables 1, 2 and 4, 5 contain the information about the P t 56 
(see formula (3) of [1]) spectrum that in reality corresponds approximately, up to the fixed 
value of (kr) = 0.44 GeV/c, to the p t ISR spectrum. 

From the comparison of Table 1 with Table 4 one can conclude that the most populated 
part of P t 56 (or P t T R ) spectrum reduces practically by twice with restricting Ptcur - So for 
A(j> max = 15° we see that it drops from < P t 56 < 50 GeV/c for P t cm = 20 GeV/c 
to a more narrow interval of < P t 56 < 20 GeV/c for the Ptcur = 5GeV/c case. 
At higher P t 7 interval for the same value A(f> max — 15° about the same factor 2 of Pt56 
spectrum reduction (from an interval < P t 56 < 50 GeV/c for Ptcur = 20 GeV/c to a 
< P t 56 < 20 GeV/c for P t cuT = 5 GeV/c) can be found from the comparison of Table 
2 with Table 5. 

From the same tables one can also see that P t 56 spectrum becomes slightly wider 
for higher values of P t 7 : the dominant part of this spectrum in the interval 40 < Pt 7 < 
50 GeV/c spreads in the case of Ptcur = 5 GeV/c cut mainly within the interval < 
Pt56 < 10 GeV/c with the maximum peak at < P(56 < 5 GeV/c. At higher PC interval 
200 < Pt 7 < 240 GeV/c, Pt56 spectrum spreads out in a wider interval of < Pt56 < 
15 GeV/c having the peak of its maximum at 5 < P(56 < 10 GeV/c. 

Thus, we can summarize that the PYTHIA generator used here predicts the values 
of p t ISR spectrum to increase while growing P t ' f , but its contribution can be reduced by 
imposing a restrictive cut on the value of p t clust (for more details see Section 3 and the 
following papers [2, 3]). 

2 The enlarging of Pt 1 values produces the same effect as it is seen from Tables 2 and 4 and will be demonstrated 
in more detail in our following papers [2, 3], 

3 We thank E. Pilon and J. Ph. Jouliet for the information about new FNAL data on this subject and on the impor- 
tance of NLO corrections and soft physics effects. 
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Table 1: Number of events dependence on Pt56 and A</> for 
40 < Pt~i < 50 GeV/c and Pt^T = 20 GeV/c for L^^/h" 1 



P t 56 
(GeV/c) 




180° 


15° 


10° 


5° 


0-5 


1 llo / /z 




lUUooz / 




5-10 








O 1ZjU4 


10-15 


1331589 


1122473 


771060 


380122 


15-20 


992374 


568279 


365329 


179767 


20-25 


725537 


282135 


183406 


91113 


25-30 


559350 


169186 


112308 


58395 


30-40 


911942 


265961 


178048 


89867 


40-50 


388950 


94112 


62068 


31000 


50-100 


91248 


19442 


12973 


6234 


100-300 


34 











300 - 500 














0-500 


7750799 


5135671 


4095348 


2498507 



Table 2: Number of events dependence on Pt56 and A<j> for 
200 < PV < 240 GeV/c and Ptcur = 20 GeV/c for L i „ t =3/b- 1 



P t 56 
(GeV/c) 




180° 


15° 


10° 


5° 


0-5 


1429 


1429 


1427 


1380 


5-10 


3266 


3266 


3264 


3150 


10-15 


3205 


3205 


3200 


3069 


15-20 


2827 


2827 


2819 


2618 


20-25 


2409 


2408 


2393 


1918 


25-30 


2006 


2006 


1982 


1300 


30-40 


2608 


2605 


2533 


1411 


40-50 


1237 


1230 


1067 


586 


50-100 


1066 


1018 


842 


536 


100-300 


313 


307 


293 


221 


300-500 














0-500 


20366 


20301 


19820 


16189 



Table 3: Number of events dependence on A0 maa; and on PC for Li nt = 3 fb 1 . 
P t gJ$ = 20 GeV/c (summary) 



{GeV/c) 




L_ 180 ° 


15° 


10° 


5° 


40-50 


7750799 


5135671 


4095348 


2498507 


100-120 


323766 


297323 


258691 


176308 


200 - 240 


20366 


20301 


19820 


16189 


300-360 


3638 


3638 


3627 


3323 
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Table 4: Number of events dependence on Pt56 and A(f> for 
40 < PV < 50 GeV/c and Pt%y% = 5 GeV/c for L^^S/fe" 1 



P t 56 




A<f> 


max 




(GeV/c) 


180° 


15° 


10° 


5° 


0-5 


331522 


331321 


329876 


295759 


5-10 


319153 


318581 


299960 


187089 


10-15 


88603 


82586 


60537 


32335 


15-20 


21244 


15327 


11663 


6924 


20-25 


8101 


5681 


4639 


2992 


25-30 


4739 


3395 


2823 


1949 


30-40 


3495 


2790 


2555 


1714 


40-50 


1647 


1277 


1042 


471 


50 - 100 


101 


67 


67 


67 


100 - 500 














0-500 


778606 


761026 


713161 


529299 



Table 5: Number of events dependence on Pt 56 and Acf> for 
200 < PtT < 240 GeV/c and PtcuT = 5 GeV / c for Uut^fh' 1 



P t 56 
{GeV/c) 


max 


180° 


15° 


10° 


5° 


0-5 


369 


369 


369 


369 


5- 10 


563 


563 


563 


562 


10- 15 


217 


217 


217 


217 


15-20 


56 


56 


56 


56 


20-25 


20 


20 


20 


18 


25 - 30 


9 


9 


9 


7 


30-40 


7 


7 


7 


6 


40-50 


6 


6 


6 


5 


50 - 100 


10 


10 


10 


10 


100 - 300 


8 


8 


8 


8 


300 - 500 














0-500 


1264 


1264 


1264 


1257 



Table 6: Number of events dependence on A<f> m ax and on Pt 7 for Li nt = 3 fb 1 . 
P t gJ$ = 5 GeV/c (summary) 



Pt 1 

(GeV/c) 




180° 


15° 10° 


5° 


40-50 


778606 


761026 


713161 


529299 


100-120 


22170 


22143 


22038 


20786 


200 - 240 


1264 


1264 


1264 


1257 


300 - 360 


212 


212 


212 


212 



4 



Since the last lines in Tables 1, 2 and 4, 5 contain an important information on A<f> 
dependence of the total number of events, we supply these tables with the summarizing Tables 
3 and 6. They include more intervals of PP and contain analogous numbers of events that 
can be collected in different Acf> intervals for two different p t clust cuts at L int = 3 (one 
month of LHC continuous running at low luminosity). 

2.2 PP and rp dependence of rates. 

Since there are two main objects for experimental regis- 
tration during the calibration procedure, namely a pho- 
ton and a jet, we shall present here the number of events 
predicted by PYTHIA simulation with cuts defined by 
(1) for different intervals of PP and rp '. The lines of 
Table 7 correspond to P t 7 -intervals while the columns 
— to rp -intervals . The last column of this table con- 
tains the total number of events in the whole ECAL (at 
Lint — 3 fb~ x ) rp -region < |?7 7 | < 2.61 for a given 
P t 7 -interval. We see from here that the events number 
decreases fastly with the Pp growth (by more than 50% 
for each following interval) while for the fixed Pp - in- 
terval there are no big changes with the variation of rp . 
Since L int = 3 fb^ 1 corresponds approximately to one 
month of LHC continuous running, we conclude that 
these rates demonstrate that, in principle, there would 



Table 7: Rates for L int = 3 fb' 1 for different intervals of PP and rf< (PtcuT = 5 GeV/c and A<j> < 15°). 



P 7 

(GeV/c) 


rp intervals 


all rp 


0.0-0.4 


0.4-0.7 


0.7-1.1 


1.1-1.5 


1.5-1.9 


1.9-2.2 


2.2-2.6 


0.0-2.6 


40-50 


102656 


107148 


100668 


103903 


103499 


116674 


126546 


761027 


50-60 


43905 


41729 


41074 


45085 


42974 


47640 


50310 


312697 


60-70 


18153 


18326 


19190 


20435 


20816 


19432 


23650 


140005 


70-80 


9848 


10211 


9963 


10166 


9951 


11397 


10447 


71984 


80-90 


5287 


5921 


5104 


5823 


5385 


6067 


5923 


39509 


90-100 


2899 


3033 


3033 


3326 


3119 


3265 


3558 


22234 


100-120 


2908 


3091 


2995 


3305 


3133 


3282 


3429 


22143 


120-140 


1336 


1359 


1189 


1346 


1326 


1499 


1471 


9525 


140-160 


624 


643 


626 


674 


706 


614 


668 


4555 


160-200 


561 


469 


557 


555 


519 


555 


557 


3774 


200 - 240 


187 


176 


186 


192 


187 


185 


151 


1264 


240 - 300 


103 


98 


98 


98 


100 


92 


74 


665 


300 - 360 


34 


34 


33 


32 


31 


27 


20 


212 


40 - 360 


188517 


192274 


184734 


194957 


191761 


210742 


226819 


1389484 



be quite enough number of events for the jet mass scale setting and successful calibration of 
HCAL. Fig. 1 gives a compact illustration of rates for four P t 7 -intervals, considered here, 

completely contained in these regions, while column HE+HF gives the number of the events 
where the jet covers both the HE and HF regions. From these tables we can see what number 



10" r 



40<P,'< 50 GeV 



100<P/<120 GeV 



200<P/<240 GeV 



300<P,'<360 GeV 
O 0.5 1 1.5 2 2.5 

Fig. 1 : ^-dependence of rates for 
different P t 7 -intervals 
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and it shows a very weak dependence on rf 1 (or better to speak about their independence on 

2.3 Estimation of "7 + Jet" events rates for HB, HE and HF regions. 

Since the jet is a wide spread object, the rf et dependence of rates for different P t 7 -intervals 
will be presented in a different way than in Table 7. Namely, Tables 8-1 1 include the rates 
of events for different rf et intervals, covered by the barrel, endcap and forward (HB, HE 
and HF) parts of HCAL and for different p t Jet intervals (for a case of Lj„ t = 3 fb^ 1 ). The 
selection cuts are taken as those of Section 3.2 of [1] specified by the following values of 
cuts parameters: 



Table 8: Selection 1. AP t jet /P t jet = 0.00 



1 p * 7 


HB 


HB+HE 


HE 


HE+HF 


HF 


40-50 


260259 


211356 


141759 


102299 


45354 


50-60 


108827 


89126 


55975 


41553 


17216 


60-70 


49585 


40076 


25172 


18153 


7019 


70-80 


25506 


20897 


12881 


9679 


3021 


80-90 


14083 


11720 


7529 


4873 


1304 


90 - 100 


7261 


7054 


4142 


2924 


853 


100- 120 


7703 


6913 


4013 


2926 


588 


120- 140 


3372 


2977 


1805 


1164 


207 


140- 160 


1650 


1481 


865 


509 


50 


160-200 


1493 


1137 


708 


396 


40 


200 - 240 


503 


406 


242 


107 


6 


240-300 


287 


215 


122 


40 


1 


300-360 


96 


73 


35 


8 





40 - 360 


480538 


393378 


255266 


184642 


75660 



Table 9: Selection 1. AP t jet /P t jet < 0.10 



Pt 1 


HB 


HB+HE 


HE 


HE+HF 


HF 


40-50 


341043 


55160 


263629 


26653 


74534 


50-60 


144955 


20396 


108765 


9300 


29281 


60-70 


65525 


8541 


49412 


3907 


12621 


70-80 


34155 


4093 


25918 


1957 


5860 


80-90 


19224 


1961 


14741 


804 


2778 


90-100 


10258 


1304 


8394 


536 


1742 


100-120 


10859 


1043 


8357 


545 


1338 


120-140 


4618 


509 


3675 


178 


546 


140-160 


2325 


222 


1751 


90 


168 


160 - 200 


1971 


147 


1458 


52 


147 


200 - 240 


685 


61 


472 


20 


26 


240 - 300 


383 


32 


234 


7 


9 


300 - 360 


129 


10 


72 


1 





40 - 360 


636418 


93480 


486788 


44052 


129050 
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P t %# T = 5 GeV/c; e^ T = 7%; < 15°; P t %$ = 5 GeV/c. (2) 

No any restrictions on other parameters are used. The first columns of these tables give the 
number of events with jets (found by LUCELL jetfinding algorithm of PYTHIA), all particles 
of which are comprised completely (100%) in the barrel part (HB) and there is a 0% sharing 
of P t jet (AP t jet ) between HB and neighbouring HE part of HCAL, i.e. AP 4 jet = 0. Second 
columns of these tables contain a number of the events in which P t of the jet is shared between 
HB and HE regions. The same sequence of restriction conditions takes place in the following 
columns. Thus, HE and HF columns include a number of events with the jets 



Table 10: Selection 2. AP t jet / P t jet = 0.00 



Pt 1 


HB 


HB+HE 


HE 


HE+HF 


HF 


40-50 


46972 


32954 


26114 


16208 


10041 


50-60 


23717 


18911 


13448 


8367 


5047 


60-70 


14384 


9751 


7469 


4703 


2386 


70-80 


8546 


6733 


4627 


2960 


1206 


80-90 


5653 


4386 


3107 


1925 


573 


90 - 100 


3326 


3119 


1900 


1377 


390 


100-120 


4157 


3435 


2271 


1467 


324 


120 - 140 


2183 


1786 


1185 


710 


134 


140-160 


1175 


1005 


635 


362 


31 


160 - 200 


1179 


905 


565 


314 


25 


200 - 240 


442 


353 


212 


97 


5 


240 - 300 


273 


200 


116 


37 


1 


300 - 360 


94 


71 


35 


7 





40 - 360 


112111 


83617 


61686 


38535 


20163 



Table 11: Selection 2. AP t jet /P t jet < 0.10 



j p * 7 


HB 


HB+HE 


HE 


HE+HF 


HF 


40-50 


60113 


7986 


45388 


3909 


14894 


50-60 


31495 


3631 


25134 


1971 


7259 


60-70 


18326 


2248 


13139 


968 


4011 


70-80 


11385 


1243 


8741 


573 


2132 


80-90 


7614 


633 


5957 


292 


1145 


90 - 100 


4544 


536 


3886 


280 


865 


100- 120 


5771 


481 


4434 


278 


689 


120- 140 


2909 


272 


2370 


94 


352 


140- 160 


1648 


138 


1246 


65 


111 


160-200 


1560 


113 


1162 


38 


115 


200 - 240 


600 


53 


416 


17 


23 


240-300 


362 


30 


220 


6 


8 


300-360 


126 


10 


71 


1 





40 - 360 


146468 


17374 


112177 


8492 


31603 
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of the events can, in principle, be suitable for the most precise calibration procedure, carried 
out separately for HB, HE and HF parts of HCAL in different intervals of Pt jet . 

Less restrictive conditions, when up to 10% of the jet Pt are allowed to be shared 
between HB, HE and HF parts of HCAL, are given in Tables 9 and 11. Tables 8 and 9 cor- 
respond to the case of Selection 1 (see [1]). Tables 10 and 11 contain a number of events 
collected with addition of Selection 2 restriction, i.e. they include only the events with "iso- 
lated jets" (whose definition is given by (26) in the same Section 3.2 of [1]). 

From Table 8, that corresponds to the most restrictive selection and gives the number 
of events most suitable for HCAL calibration, we see from the last summarizing lines that for 
the interval 40 < PV < 360 GeV/c PYTHIA predicts around a half of million of events for 
HB and a quarter of million of events for HE per one month of continuous data taking at low 
LHC luminosity, while for HF it is expected to be less than 80 000 events per month. 

One should keep in mind that the last columns in Tables 8-11 can not be taken as the 
final result here as we have not defined the meaning of sharing Pt 3et between the HF regions 
and the region with \q\ > 5, i.e. close to a "beam-pipe" region. A more accurate estimation 
can be done here by finding events with jets in a wider region than the HF volume restricted 
by 3 < \rj HF | < 5 and by calculation of a number of those events whose jets are contained 
in HF completely with all particles belonging to them. 

3. STUDY OF FEATURES OF "7 + Jet" EVENTS IN THE HCAL BARREL RE- 
GION. 

In this section we shall study a specific sample of events that may be most suitable for HB 
calibration. It is a sample of events in which jets are completely (100%) contained in HB 
region, i.e. having 0% sharing of Pt jet (at PYTHIA level of simulation) with HE. Below 
we shall call them "HB-events". A particular set of these events for p t clust = 5 GeV/c 
is presented in the first column (HB) of Table 8. Here we shall use two different jetfinders, 
namely, LUCELL from PYTHIA and UA1 (the last one is taken [] from CMSJET [6]) for 
an equal foot determination of clusters and jets. The distributions of p t clust f or generated 
events found by these two different jetfinders in two P t 7 intervals, 40 < Pt' < 50GeV/c and 
300 < Pt 1 < 360 GeV/c, are shown in Fig. 2 for P t $$ = 30 GeV/c. 

3.1 Influence of the Ptcur parameter on the photon and jet P t balance and initial 
state radiation suppression. 

Here we shall study correlation of p t clust with p t ISR proportional to P t 56 for the fixed (k T ) 
(The influence of the (kx) variation on the P t 7 and P t Jet balance is discussed in [4]). 

The banks of 1-jet "7 + Jet" events gained from the already mentioned (in Section 
2.1) results of PYTHIA generation of 10 7 signal "7+ Jet" events in each of four P t 7 intervals 
(40 - 50, 100 - 120, 200 - 240, 300 - 360 GeV/c) will be used here. The observables defined 
in Sections 3.1 and 3.2 of paper [1] will be restricted here by Selection 1 cuts (17) - (24) of 
Section 3.2 of [1] and the cut parameters defined by (1). The luminosity was chosen to be 
L = 10 33 cm 2 sec~ 1 . 

4 In order to select events with p t clust value starting from 5 GeV/c we have changed the Pt precluster threshold 
in the UA1 algorithm from 10 GeV/c, taken as default value, to 1.5 GeV/c . We also increase the cone radius in 
the UA1 algorithm from 0.5 to 0.7. 
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We have chosen two of these intervals with the extreme values of Pp to illustrate the 
influence of the PtouT parameter on the distributions of physical variables. The results of 




Fig. 2: Pt c! " s '-distribution in "7 + Jet" events from two Pt 1 intervals: a) 40 < Pt 1 < 50 GeV/c and 
b) 300 < Pf < 360 GeV/c with the same cut Ptcur = 30 GeF/c. 

Ptcuf variation are shown in Fig. 3 for 40 < P t 7 < 50 GeV/c and in Fig. 4 for 300 < 
Pt 1 < 360 GeV/c. In these figures, in addition to three variables P t 56, P t v>5 , P t out , 
already explained in Sections 2.2, 3.1 and 3.2 of [1], we present the distributions for other 
two variables, P t (0+f] > 5) and (1 — cosAifi), which define the right-hand side of equation 
(29) of [1]. The distribution for the back-to-back A0 angle (see (23) of [1]) which defines 
the second variable, is also presented in Figs. 3, 4. 

The disbalance variable Pf56 (defined by formula (3) of [1]) and both components of 
another disbalance measure (PJ - P/ et )/P{i (defined by formula (29) of [1]) as a sum of 
(1 — cosAcj)) and P t (0 + r] > 5) , as well as two others, p t out and A</>, show a tendency, as 
it is seen from Figs. 3 and 4, to become smaller with the restriction of the upper limit on the 
p^ciust va j ue i t means that the calibration precision may increase with decreasing PtcuT' 
which justifies the intuitive choice of this new variable in Section 3 of [1]. The origin of this 
improvement becomes clear from the P t 56 density plot, which demonstrates suppression of 
P t 56 (or P t ISR ) with implying a more restrictive cut on p t clust . 

The comparison of Fig. 3 (for 40 < P t 7 < 50 GeV/c) and Fig. 4 (for 300 < PC < 
360 GeV/c) shows that a degree of back-to-backness of the photon and jet P t vectors in 
the 0-plane increases with increasing Pf . At the same time p t out and p t ISR distributions 
become wider, while the Pt v>5 distribution practically does not depend on Pf and p t clust _ 

It should be mentioned that the results presented in Figs. 3 and 4 were obtained with 
the LUCELL jetfinder of PYTHIA 0. 



5 The results obtained with both jetfinders and Pt 1 and p t Jet balance will be discussed in [2, 3] in more detail 
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Fig. 3: LUCELL algorithm, Acf> < 15°, 40 < Pt< < 50 GeV/c. Selection 1. 
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Fig. 4: LUCELL algorithm, Acp < 15°, 300 < Pt'< < 360 GeV/c. Selection 1. 
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3.2 Jetfinders and P t structure of jets in ri - <f> space. 

In order to understand well the calibration procedure of "7 + Jet" events, it is useful to keep 
control over some principal characteristics of jets and of the P t activity in the space around 
them. The left-hand columns of Figs. 5 and 6 represent jet radius W et (r), (j>) distributions in 
these HB-events for the intervals: 40 < P t 7 < 50 GeV/c and 300 < Pt 1 < 360 GeV/c, 
respectively. We have chosen the jet radius counted from the initiator cell (ic) to be restricted 
by Rlf = 0.7 for the LUCELL and UA1 jetfinders. 

From the left-hand side plots in Figs. 5 and 6 we see that both jetfinders give close 
enough results for distribution. In these plots the radius in r\ — <fi space is counted from 
the gravity center (gc). The detailed information about averaged values of the jet radius for 
four Pt 7 -intervals will be presented in the tables of Appendices \-A of [3] []. 

Let us now consider the question, how the transverse momentum is distributed inside 
a jet. Let us divide the jet radius W et {j] : <j>) = R into a set of AR bins and calculate the 
vector sums of cells P t in each ARu n ring. Normalized to P t Jet , the modulus of this vector 
sum, denoted by P t bm , would give the value that tells us what portion of P t Jet is contained 
in the ring of size ARbi n . Its variation with the distance R counted from the center of gravity 
of the jet is shown in the right-hand columns of Figs. 5 and 6. 

From these figures we can conclude that the LUCELL and UA1 jetfinders lead to 
equally P f -densed central part inside the jet. 

3.3 P t distribution inside and outside a jet. 

Now let us see how the volume outside the jet, (i.e. calorimeter cells outside the jet cone), 
may be populated by P t in these HB "7 + Jet" events. For this purpose we calculate a vector 

-* sum 

sum P t of individual transverse momenta of Arj x A(p cells, included by a jetfinder into 
a jet as well as of cells in a larger volume that surrounds a jet. In the last case this procedure 
in some sense can be viewed as some straightforward enlarging of the jet radius in the r\ — <fi 
space. 

The figures that represent the ratio P t sum / P^ , as a function of the distance R(r], (f) 
counted from the jet gravity center towards its boundary and further into space outside the jet 
are shown in the left-hand columns of Figs. 7 and 8 for a case when all jet particles are kept 
in the jet, while the case of the magnetic field effect account (i.e. when the contribution from 
charged particles with Pt < 0.8 GeV/c is removed ^ from the total jet Pt in a case of HB 
events) is shown in the left-hand columns of Figs. 9 and 10. 

From these figures we see that the jet surrounding space is found to be far from being 
an empty one in the case of "7 + Jet" events considered here. We also see that an average 
value of the total p t sum increases with increasing the volume around the jet and it exceeds 
Pt 1 at R = 0.7 — 0.8 when all particles are included into the jet (see Figs. 7 and 8), while in 
the case of rejected charged particles with P t ch < 0.8 GeV/c only about 95 - 97% of P t Jet 
are collected at R = 0.7 - 0.8. 

From the right-hand column of Fig. 7 we see that when all particles are included into 
the jet, the following disbalance measure: 



6 From there one can see a weak dependence of the jet radius on Pt Jet (m Pt 1 ) for both UA1 and LUCELL 
algorithms 
7 See [5] 
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Fig. 5: LUCELL and UA1 algorithms, A</> < 15°, 40 < Pi 7 < 50 GeV/c. Selection 1. 
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Fig. 6: LUCELL and UA1 algorithms, A0 < 15°, 300 < Pi 7 < 360 GeV/c. Selection 1. 
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Fig. 7: LUCELL and UA1 algorithms, A0 < 15°, 40 < Pt< < 50 GeV/c 
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Fig. 8: LUCELL and UA1 algorithms, A</> < 15°, 300 < Ft 1 < 360 GeV/c 



16 




17 




Fig. 10: LUCELL and UA1 algorithms, Acfi < 15°, 300 < Pt 7 < 360 GeV/c, P t 3 ^ > 0.8 GeV/c. 
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p 7+s-um 



-+ 7 -* sum 

P t + P t 



(3) 



achieves its minimum at P « 0.7 — 0.8 for both LUCELL and UA1 algorithms. From 
Figs. 9 and 10 we see that with account of the magnetic field effect on charged particles with 
P t ch < 0.8 GeV/c the minimum of pf+ sum i s achieved at larger values of R « 1.0 again 
for both LUCELL and UA1 jetfinders. 

The value of p t "i+ sum continues to grow rapidly with increasing R after the point 
R=\ for 40 < Pt 1 < 50 GeV/c (see Figs. 7, 9), while for higher P t 7 (see Figs. 8, 10 for 
the 300 < Pt 1 < 360 GeV/c interval) the ratio p t sum /p t ~< and the disbalance measure 
p^i+sum mcrease more slowly with increasing R after the point R=0.7. This means that at 
higher P t 7 (or Pt' Jet ) the topology of "7 + Jet" events becomes more pronounced and we 
get a cleaner picture of an "isolated" jet. This feature clarifies the motivation of introduction 
of two criteria "Selection 2", (see "Point 8") and of "Selection 3" (see "Point 9") in Section 
3.2 of [1] for selection of events with "isolated jets". 



4. SUMMARY 

From the study in Section 2.1 we see that the PYTHIA simulation predicts that most "7 + 

-* 7 -* Jet 

Jet" events at LHC energies may have the vectors P t and P t being back-to-back within 
Acj> < 15° (about 66% for 40 < P t 7 < 50 GeV/c and greater than 99% for P t 7 > 
200 GeV/c). At the same time, as it is seen from Table 1, a substantial part of the P t 7 (« 
Pt 56) spectrum spreads out in the interval < Pt 56 < 20 GeV/c, having the peak inside 
the 5 < Pt 56 < 10 GeV/c interval and noticeable tails (about of 5 — 10% of number of 
events) extending to P(56 = 40 GeV/c not only for A<p < 15° but also for the smaller 
A</> < 10° and A<j> < 5° intervals. From here we can conclude that the use of only one 
traditional A<f> cut does not help much to reduce the p t ISR contribution at LHC energies. It 
is most obvious from Table 2, where we find that at high P t 7 (200 < P t 7 < 240 GeV/c) 
more that 99% of events belong to the A<j) < 15° interval. Most of these events lay within 
the < P(56 < 50 GeV/c interval. Even for Ac/) < 5°, a considerable part of the spectrum 
belongs to the interval < P*56 < 40 GeV/c. Thus, a decrease in A<j> has, in fact, no big 
influence on the p t ISR spectrum and it cannot be useful to reduce the number of events with 
essential amount of initial state radiation at LHC energies. 

In Section 3 the efficiency of Ptcur f° r ISR suppression was demonstrated. From 
Figs. 3 and 4 we see that the distribution for P(56 becomes narrower with decreasing Ptcur- 
Analogous behavior is shown by the P t out (as p t clust i s a part of p t ° ut ) spectra. These 
figures serve as an illustration for a more detailed study of Ptcur influence on P t 7 and 
P t Jet balance, which will be presented in our next papers [2-4]. 

A strict limitation of the Ptcur parameter, introduced in Section 3.2 of [1], (by w 5 — 
10 GeV/c) improves this tendency. Simultaneously, due to this limitation, one can noticeably 
reduce radiation in the initial state (compare Tables 1, 4 and 2, 5), which leads to decreasing 
P t 7 and P t Jet disbalance. 

From Section 2.3 (Table 8) we see that even with the restrictive cuts mentioned there 
one can expect around a half million of events for HB, a quarter million of events for HE, and 
less than 80 000 events for the HF part and Li nt = 3 fb^ 1 (i.e. per month of continuous data 
taking at low LHC luminosity) for interval 40 < Pi 7 < 360 GeV/ c for the sample of events 
that corresponds to 0% P t sharing between HCAL parts. Despite those cuts the number of 
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events for the HCAL part above, in principle, can be quite sufficient for successful jet energy 
scale setting and HCAL calibration. 
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